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Abstract 

It is demonstrated that radiative corrections increase tunneling probability 
of a charged particle. 
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Famous baron von Miinchhausen saved himself from a swamp pulling his hairs by the 
hands of his own According to classical physics, such a feat seems to be impossible. 
However, we live in a quantum world. In a tunneling of a charged particle, the head of 
the particle wave function can send a photon to the tail which absorbs this photon and 
penetrates the barrier with enhanced probability. Obviously, such a photon feedback should 
work in the two-body tunneling where the first particle, while continuing to be accelerated 
by the potential after the tunneling, can emit a (virtual) photon that increase energy of the 
second particle and its tunneling probability. The Miinchhausen mechanism may be helpful 
in the tunneling of a composite system. It is related to phonon assisted tunneling but does 
not require any special device being always provided by the interaction of a charged particle 
with the radiation field. 

The interaction of a tunneling object with other degrees of freedom of the system and 
the influence of this interaction on the tunneling probability for a long time was a topic 
of intensive studies initiated by Caldeira and Leggett 0. Their general conclusion, in 
agreement with intuitive arguments, was that any friction-type interaction suppresses the 
tunneling. At the same time, it was realized that such an interaction leads to distortions of 
the barrier which can be helpful in endorsing the tunneling. The simplest effect is associated 
with the zero-point vibrations of the source responsible for the existence of the barrier. 
This is important for the probabilities of subbarrier nuclear reactions as pointed out by 
Esbensen || . In the last decade, many experimental and theoretical efforts were devoted to 
the understanding of related aspects of subbarrier reactions, see the recent review [|J and 
references therein. Below we discuss the interaction of a charged tunneling object with the 
electromagnetic field which always accompanies motion of the object. 

Formally speaking, we are looking for the effects of radiative corrections on the single- 
particle tunneling. These effects can be described by the Schrodinger equation with the 
self-energy operator: 

HV(t) + f E(r, r'; £)*(r')dV = E^(r) (1) 
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where H is the unperturbed particle hamiltonian, which includes a barrier potential, and 
X = M — iT /2 is the complex nonlocal and energy- dependent operator determined by 
the coupling to virtual photons and a possibility of a real photon emission. The "photon 
hand" here connects two points r and r' of the same wave function. In the one-photon 
approximation the self-energy due to the interaction with the transverse radiation field can 
be written as 

Here we sum over unperturbed stationary states |n); f and p are the position and momentum 
operators, respectively; the photons are characterized by the momentum k, frequency and 
polarization A; the polarization vectors are perpendicular to k so that the momentum 
operators commute with the exponents. The normalization factors are included into oc 
Lu k . The relativistic generalization of (|2|) is straightforward. 

The hermitian part M of the self-energy operator is given by the principal value integral 
over photon frequencies in (Q). The expectation value of M is responsible for the Lamb 
shift of bound energy levels. It contains also the mass renormalization for a free particle 
which should be subtracted. Our problem is different from the energy shift calculation for 
bound states since we are interested in the change of the wave function of the tunneling 
particle. However we can use some features of the conventional approach. As well known 
from the Lamb shift calculations, one can use different approximations in the two regions 
of integration over the photon frequency to. In the nonrelativistic low- frequency region, 
uj < [3m, where the parameter (3 < 1 is chosen in such a way that typical excitation energies 
of a particle in the well SE are smaller than (3m (in the hydrogen Lamb shift problem a 
fine structure constant a can play the role of the borderline scale parameter), it is possible 
to neglect the exponential factors in (fj). The high-frequency contribution to M, where the 
potential can be considered as a perturbation to free motion, has been calculated, e. g. in 
Ref. [Q. The two contributions match smoothly at u — (3m. 

It is easy to estimate the mass operator M with logarithmic accuracy. After summation 



over polarizations and standard regularization || , the low frequency part of the operator M 
can be written as 

where Ze is the particle charge, and m is the mass of the particle (reduced mass in the 
alpha-decay case). We use the units % = c = 1. Substituting the logarithm arising from the 
frequency integration by its average value L = \n.(/3m/u) min ), we can use the closure relations 
and obtain a simple expression 

M(E) = f^iP(ff " E)i> (4) 

= £% L {™ + K A (5) 

The mean value of the term with the anticommutator [...,...] + in eq. (|5]) is equal to zero 
since (H — E)^o = where \I/o is the unperturbed wave function. A correction to the 
wave function due to this term can be calculated by using perturbation theory and the 
unperturbed Schrodinger equation, 

**=^Lp7-(0|E/|0)]*o. (6) 
Sum 

This correction is not essential since it does not influence the exponent in the tunneling 
amplitude. 

Combining the remaining term in eq. (|5|) with the high-frequency contribution which 
contains L = ln(m/ /3m), see Ref. 0, the result can be presented as an effective local operator 
proportional to the Laplacian V 2 U(r), 

M(r, r'; E) ~ V 2 U(r)5(r - r ')—^ln — = 5U(r)5(r - r'). (7) 

3nm 2 U 

Here we used the barrier height U as a lower cut-off u) min of the integration over frequencies 
(below we give a semiclassical estimate which leads also to a more accurate evaluation of the 
logarithmic factor). For the tunneling of an extended object, the mass m in the argument of 
the logarithm should be replaced by the inverse size of the particle l/r which comes from 



the upper frequency cut-off given in this case by the charge formfactor. The obtained result 
is physically equivalent to the averaging over the position fluctuations due to the coupling 
to virtual photons. Thus, in the logarithmic approximation the mass operator is reduced to 
a local correction 5U(r) to the potential U(r). 

The Laplacian of the potential energy V 2 U(r) near the maximum of the barrier is neg- 
ative (correspondingly, near the bottom of the potential well it is positive). Therefore, we 
obtained the negative correction SU (r) to the potential barrier which leads to a conclusion 
that jiggling of the photon increases the tunneling amplitude of the particle. The numerical 
value of the correction to the potential is small, ~ 1 keV, for the alpha-decay. However, in 
some cases it may be noticeable due to the exponential dependence on the height of the bar- 
rier (recall the notorious cold fusion problem). Also, there exist theories like QCD where the 
radiation corrections are not small. In many-body systems one can use collective modes, as 
phonons, to transfer energy. This can influence electron tunneling through quantum dots or 
insulating surfaces. Energy exchange between a tunneling particle and nuclear environment 
is known to be important in subbarrier nuclear fission and fusion ||. 

An analysis can be performed more in detail by using the semiclassical WKB approxima- 
tion for the tunneling wave functions. The semiclassical radial Green function of unperturbed 
motion under the barrier can be written in terms of the classical momentum in the forbidden 
region, p(r; E) = [2m(U(r) — E)] 1 ^ 2 , at a given energy E as 




where Q(x) is the step-function. The full three-dimensional Green function G(E) = 
J2 n \ n ){E — E n )~ l (n\ contains also angular harmonics which could be separated in a routine 
way accounting for the fact that in the long wavelength approximation for the s-wave solu- 
tion the intermediate states are p-waves. Indeed, the operator of electric dipole radiation p 
converts an initial s-wave $ in eq. ([I]) into an intermediate p-wave state \n). Therefore, it 
is sufficient to keep the p-wave part of the radial Green function and to use closure in the 
sum over angular harmonics. 
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The kernel of the integral term in the Schrodinger equation (1) contains 

K(r, r'; E) = J duo G(r, r'; E-u). (9) 

The integrand consists of terms falling exponentially as \r — r'\ increases. The potential U(r) 
is assumed to be a smooth function. Therefore, we can put p(r') ~ p(r). Now it is easy to 
perform the integration over uj in eq.@ which leads to 

K(r,r';E) = - 1 — { e ~P™n\r-r'\ _ e -p ma ,\r-r'\ j ^q) 

where p m in = [2m(U p (r) — E)] 1 / 2 , p ma x = [2/?] 1//2 m, and U p (r) is the effective p-wave radial 
potential which includes the centrifugal part. This expression has a very narrow maximum 
near r = r' with the width \r — r'\ ~ l/Pmax- This is a measure of non-locality of the 
self-energy operator M(r, r'\ E). In any nonrelativistic application the kernel can be treated 
as proportional to the delta-function. The proportionality coefficent can be found by the 
integration over r. Thus, we obtain the local behavior of the kernel, 

K(r,r';E) w -L(r)5(r -r'), (11) 

where now we determine the lower limit of the logarithm which has appeared in our previous 
derivation (0) as related to the local value of the potential, 

L(r) = l »\m=E\ (12) 

The substitution into eq. (0) gives 

su ^-^mrE\^ wy (13) 

As usual, this semiclassical expression is not valid near the turning points where U p (r) = E. 
However, a very weak logarithmic singularity does not produce any practical limitations on 
the applicabilty of eq. (p~3|). 

The conclusion of enhancement of the tunneling probability seems to contradict to the 
common sense: radiation should cause energy losses and reduce the tunneling amplitude 
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of the charged particle. However, such an argument may be valid only for the real photon 
emission. This emission is described by the antihermitian part of the self-energy operator 
which is originated from the delta-function corresponding to on-shell processes, 

r(r,r';£) = [ dcu ^(r|pe ik >) (n|p e - lkf \r')u8{E - E n - to). (14) 

Because of the energy conservation the sum here includes only states \n) with energy E n 
below E. Consider for example the tunneling from the ground s-state. A dipole transition 
transfers the particle from the s-state to a p-state. However, there are no quasidiscrete p- 
states |n) below the ground state in the potential well. Scattering p- waves can penetrate the 
potential barrier from the continuum with an exponentially small amplitude. This means 
that T(r, r'; E) is again exponentially small if one or both arguments r and r' are under the 
barrier or inside the potential well. Therefore, T(r,r';E) does not considerably influence 
the tunneling amplitude. The reason for that can be easily understood. A real radiation 
would be impossible if there were no tunneling. Whence, the radiation width must vanish 
together with the tunneling width. On the contrary, the real part of £ under the barrier 
would present even if the tunneling probability would vanish. 

To avoid misunderstanding we need to stress that the contribution to the radiation 
intensity from the barrier area and the potential well, which was, in application to the 
nuclear alpha-decay, the subject of recent experimental || and theoretical studies, still 
may be important. The radiation amplitude with E s — E p = uo contains the matrix element 

1 i 
< s \p\p >= — < s\\H, p]\p >= — < s\VU\p > . (15) 

When one moves inside the barrier from the outer turning point inwards, the resonance s- 
wave function exponentially increases while the non-resonance p-wave function exponentially 
decreases. As a result, the product V's( r )' ? /'p( r ) does not change considerably. This means 
that the contribution to the real radiation from the inner area may be comparable to that 
from the area outside the barrier. The gradient VU changes its sign near the maximum of the 
potential which implies a destructive interference between the radiation from the different 
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areas (since \s) is the non-oscillating ground state wave function, the product ip s (r)tp p (r) 
does not change sign inside the barrier). 

This work was supported by the NSF grant 96-05207. V.F. gratefully acknowledges the 
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